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SUPPORTING  SURFACE  IN  NONSTATIONARY  FLOW  NEAR  SCREEN 


B.  N.  Belousov,  A.  N.  Lukashenko,  A.  N.  Panchenkov 

At  the  present  time  there  exist  a number  of  Important  studies 
on  the  theory  of  a supporting  surface  in  a nonstationary  flow  [1, 

2,  3]  with  different  approximate  methods  of  calculating  the  wing. 

As  we  know,  the  problem  of  a thin  supporting  surface  of  arbitrary 
length  is  reduced  to  two-dimensional  integral  equations  for  which 
there  are  no  closed  solutions. 

One  of  the  most  Important  aspects  is  that  of  studying  two- 
dimensional  Integral  equations  in  one-dimensional  approximations 
and  obtaining  on  this  basis  final  results  for  the  hydroaerodynamic 
characteristics  of  a wing  in  a nonstationary  flow. 

Whereas  for  a wing  in  a stationary  flow  Prandtl's  theory  of 
the  supporting  line  defines  the  Prandtl  equation  in  terms  of  the 
one-dimensional  approximation,  for  a wing  in  a nonstationary  flow 
the  complexity  of  the  physical  phenomenon  does  not  provide  a unified 
approximation  method.  Thus,  there  are  a large  number  of  approaches 
(up  to  20,  as  noted  by  R.  A.  Welspllnghoff , H.  Ashley,  and  R.  A. 
Halfman  [1])  which  give  different  one-dimensional  equations. 

Developed  in  the  works  of  one  of  the  authors  [6,  7,  8,  9]  is  a 
method  of  Integral  operators  for  the  theory  of  a lifting  surface  in 
a nonstationary  flow. 
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The  general  solution  to  the  problem  is  represented  in  the  form 
of  three  components.  Two  components  are  found  from  equations  which 
correspond  in  form  to  the  equations  of  the  theory  of  a wing  in  a 
nonstationary  flow. 

Also  examined  in  [8]  is  a method  of  constructing  the  Prandtl 
equation  based  on  the  example  of  a wing  in  an  unlimited  fluid  flow. 

A significant  advantage  of  this  method  is  that  the  Prandtl  equa- 
tion is  constructed  for  the  parametric  constant  of  a singular  solu- 
tion. This  has  made  it  possible  to  more  deeply  analyze  the  problem 
and  critically  examine  the  existing  theories  of  Relsner  [1],  Kusner 
[2] , and  others . 

In  the  present  article  a method  of  Integral  operators  is 
applied  to  the  problem  of  a supporting  surface  near  a screen. 
Obtained  are  two-dimensional  integral  equations  and  a one- 
dimensional  Intego-dlfferentlal  equation  for  the  parametric 
constant.  It  should  be  mentioned  that  the  downwash  is  associated 
only  with  the  singular  solution  corresponding  to  a singularity  in 
the  leading  edge  of  the  wing.  Transfer  of  this  result  to  the 
theory  of  stationary  motion  indicates  that  the  classical  theory  of 
a lifting  line  is  strictly  applicable  only  to  a plane  lifting 
surface . 

If  the  lifting  surface  is  not  planar,  then  in  the  expression 
for  the  coefficient  of  the  lifting  force  of  the  wing  there  appear 
terms  on  the  order  of  1/X  which  are  not  considered  by  the  classical 
theory.  From  the  physical  standpoint  they  determine  the  effect  of 
aspect  ratio  on  the  angle  of  zero  lift  of  the  wing. 

The  present  study  adheres  very  closely  to  studies  [6,  7,  8»  9], 
for  which  reason  we  have  attempted  to  very  briefly  introduce  the 
general  results  from  these  studies  in  the  necessary  cases  and 
discuss  new  results  in  more  detail. 
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1.  Let  us  assume  that  a thin  lifting  surface  S moves  at 
velocity  Vq  near  a solid  screen.  In  additional  to  the  main  forward 
motion  it  performs  harmonic  oscillations  at  frequency  w. 

Assuming  that  the  disturbances  introduced  into  the  flow  by 
the  lifting  surface  are  small,  in  coordinate  system  oxyz  moving 
at  velocity  Vq  (see  the  figure)  let  us  use  a known  method  to 
linearize  the  problem  [4,  5]. 


In  the  case  of  harmonic  oscillations  for  hydrodynamic  poten- 
tials of  velocities  cp  and  accelerations  0 and  their  derivatives 
we  have 

(1.1) 

The  index  n denotes  the  derivative  of  the  n-order  of  functions  0 
and  q>  with  respect  to  coordinate  x^^Cx^  = x;  X2  = y;  x^  = z). 

In  equations  (1.1)  and  henceforth  we  shall  only  study  the 
real  parts  of  the  proper  expressions,  although  we  shall  drop  Re 
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for  simplicity.  Relationships  (1.1)  make  it  possible  to  formulate 
the  problem  for  0(g)  and  cp(g)  and  exclude  time  from  our  study. 

By  shifting  to  dimensionless  space  assuming  that  the 
velocity  of  the  impinging  flow  equals  1 and  removing  the  bar  from 
the  top  of  corresponding  0(g)  and  tp(g)  expressions,  we  get  the 
boundary  problem  in  the  space  of  the  potential  of  the  accelera- 
tions 0: 


Ae=0;  g€Q; 
e,=F(g);  g€^y. 


0,  = 0;  ' z = 2S; 

0+  — 0_  =0.  giLi 

0-^0; 

X -►  + CO. 


(A) 


(1.2) 


(1.3) 


Problem  (A)  is  obtained  by  integral  operator  , assigned  in 

2 

space  L^(Sp)  with  values  of  C Q.  The  required  properties  of 
operator  Ay  are  given  in  [5-8]. 


In  the  studied  problem  operator  Ay  has  the  form  of 


“ S J T (/») sif  + sign F 1] dS; 

here 

ri  = K(*-£)*  + (i^-s)*  + U + C + 4A)«; 


sign  P = 


-1  - solid  screen 
+1  - free  surface. 


According  to  condition  (1.2)  we  get  the  two-dimensional 
Integral  equation  of  the  problem 


Let  us  represent  the  solution  in  the  form  of  three  components 
[6-8]: 


T (e)  = 7i(g)+  («)  + 73  («) 


where  is  the  solution  of  class  C^CSp),  related  to  the  presence 
of  a break  in  tangent  velocities  during  transfer  of  the  surface 
S J Yq  ~ solution  of  class  C^(S  ),  which  describes  inertial  move- 
ments ; y^  “ singular  solution. 

Analysis  of  equation  (1.5)  and  the  method  of  constructing 
equations  for  Yj^  is  given  in  [6], 


Let  us  write  equations  for  y-^  ^nd  y^: 


~ ^1  4-C4  g^Sj,; 
^oi-4,7i'  = —ikFt,  g e S„. 


(1.6) 

(1.7) 


Here  F2(g)  - the  function  which  describes  the  shape 


X 

of  the  surface  S;  j dx . 


Constant  C is  found  from  the  condition  of  solvability  of 
equation  (1.6)  in  the  space  of  C^(Sp)  [6]. 

By  performing  calculations  for  Nq^A^Y  we  get  [10] 


+1-H 


(x'*:  ” =2,  J J ^ IPZT)  — : 7;  [ q.v 


^ : 


^ l^-  (i-8> 


Where 


7’(P)  = ^ 


»(1)* 
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The  singular  solution  is  determined  from  equation 


Ar= 

k is  the  Strouhal  chord  number. 


(1.9) 


Equation  (1.9)  is  obtained  by  mapping  equation  (1.5)  in  the 
space  of  the  velocity  potential. 


2.  Now  let  us  examine  the  Prandtl  equation  for  a wing  and  a 
nonstationary  flow.  For  a wing  of  great  aspect  ratio  operator 
N„,A  Y is  approximated  by  operator 

U X z 


+1 

^ J T (5;  y)  [(i-rj)  + ^ 


(x-£) 


Ids. 


(2.1) 


For  a plane  flow  the  form  of  the  singular  solution  y^  is 
known,  and  for  a wing  of  great  aspect  ratio  we  use  y^  i*^  form 
of  the  solution  of  the  planar  problem,  although  with  a constant 
which  depends  on  the  y-coordinate : 


If 

JC 

tJ,  € C*  (S)  and  t*  = — ^ e 

then  we  get  [6] 

g€S„.  (2.2) 

By  examining  equations  (1.6),  (1.7),  (1.9)  and  property  (2.2), 
we  get 

*=  — C — tAsJu- 

Equation  (1.6)  can  be  written  in  the  form 


6 


NoA,^^~=F^  + 0, 
A/.  = -Sdt. 


(2.3) 


When  X ^ C = C,  while  for  greater  aspect  ratio  we  use  the 
approximation  +i 

C = — \ (2.4) 


—1 


-1 


Then,  by  Introducing  operator  p : 

P~W,.;4,7.  = f Ti(0d5 


and  taking  approximation  (2.1)  from  (2.3),  we  get  the  Prandtl 
equation  for  the  regular  solution,  which  in  form  coincides  with 
the  Prandtl  equation  of  the  stationary  problem 

' ■ w - ^ (<■■  --  “ } (2.5) 


where 


+1 


/.=]■  T.(«)d5. 

— I 

Function  Y22  found  from  equation 

= ikC(y). 


(2.6) 


Now  let  us  introduce  a number  of  assumptions  typical  for  the 
Prandtl  problem  [4],  [8]: 

N 1^*7  ^ ^ ^ 

7(P)  = 7iO)7t(P): 

^10^,7  = 271  (J^)  \{y)  (t); 

Nix  = — i.’'** 

''■At  -s  (|T(o[i;lg  + + 

■ft  » 

f 7 (E)  f r-'*’  [.— i-r.  + sign  f . 1 * dl\ 
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then 


a {y)  2)^  (y)  (0  = -C(y)-  2X  (y)  (y)  _ 

— A/uA,  |o  (y)  ]/^  + 7m]  : « € Sp. 


For  a wing  in  an  infinite  flui.  both  parts  of  equation  (2.7) 

i k X i U"  Y 

contain  factors  e . By  multiplying  equation  (2,7)  by  e~  , we 

eliminate  the  variable  x [5,  8].  This  also  occurs  in  the  general 

i kx 

case,  although  the  equation  relative  to  e cannot  be  solved  in 
explicit  form,  since  it  is  possible  to  calculate  the  integrals  in 
^lO^z''^  in  a closed  form.  In  this  connectic-  let  us  introduce 
averaged  equation  (2,7)  according  to  the  following  rule: 


For  "^222  ^'^^^10^z''^221  select  representations  in  the 

form  of 


ydft) 


72»(1/)  = 


Tw»(y) 


-Vidk;  S)2« 


(2.8) 


H)7k(y)X°(y) 


^iN,{k)C{y)  + C^(!,)N,m].  (2.9 


When  h -►  ® representations  (2.8)  and  (2.9)  are  transformed 
into  the  expressions  found  in  [91.  Functions  N^(k)  are  also 
obtained  in  [9]. 


Then , 


a{y)  2>  (y) 


+ 


1m(y) 


4-,  (*;  h)  2»  ■ 2itil/,  (*i  i)  7^^  (y) 

= — C(y)-A(jx">l,|o  (y)  ]/ 


j.|^A^,(A)C^)  + C,(y)A/,(ft)] 


(2.10) 
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+1 


with  the  Introduction  of  a new  variable  a,(y) 


+ y22i(yj 


j s)j 


dy  equation  (2.8)  is  transformed  into 


ai{y)=^  l“*+  iv,(fc)  )• 


(2.11) 


Here 


«1C(v)  + 


^ +, 


t' 

4^«  ^***  J I 


As  follows  from  the  solution  to  the  plane  problem  of  the 
nonstationary  motion  of  a wing  near  a screen  [7,  8],  function  y{S) 
can  be  written  as: 

7 (S)  = <1.0  a.  (y)  + I-tI  (S)  + (S). 

Functions  y^(S),  which  correspond  to  the  plane  problem  for  an 
unlimited  fluid,  have  been  calculated  in  [9]-  By  using  these 
results,  we  get 

^ (i')  + (*)]  (^ ->)  + (y)^.  (2.12) 

where  Laj^  is  the  integral  operator,  whose  core  was  obtained  in 
monograph  [4]. 


Now  we  must 
which  exclude 


establish  the  connection  between  the  functions 
and  <p^  from  (2.9).  By  definition  when  A -►  <» 


ax 


i'caiUf) 

~2f(yr 


J- 
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Constants  C-,(y)  take  the  form  of 


_ VT^  < « /nr: 


From  equation  (2.9)  we  find 


al 


5;^[;*.(»)Cfe,)  + C.brt».,*,](fc_,). 


By  comparing  the  two  expressions  for  a^(y),  we  get  the 
relationship  we  need 

tn  =*  ♦*;  “ to. 

Now,  combining  the  results  of  equation  (2.9),  we  write  the 
intego-differential  equation  for  a wing  of  great  aspect  ratio  near  a 
screen  in  a nonstationary  flow 

■ +1 


where 


The  core  K (9)  has  the  form  of  [4] 


10 


Where  K -*■  0,  N^Ck)  -1  and  from  (2.11)  we  get  the  Prandtl  equation 
for  a wing  near  a screen  in  a stationary  flow  Here 

Thus,  the  Prandtl  problem  for  a wing  in  a nonstationary  flow 
is  reduced  to  two  one-dimensional  Integral  equations  (2.5),  (2.13), 
where  equation  (2.5)  corresponds  to  the  equation  of  the  stationary 
theory . 

This  result,  which  is  the  direct  result  of  breaking  the 
general  solution  down  into  three  components,  is  significant  in 
that  the  authors  of  known  theory  [1,  2],  by  introducing  integral 
equations  of  the  "reduced  circulation"  type,  have  found  one  one- 
dimensional equation. 


3.  In  determining  aerodynamic  coefficients  C and  C a good 

y X 

first  approximation  gives  us  the  solutions  introduced  in  [3,  4]. 

r(i/)  = ® vT-fi. 

We  use  this  representation  for  the  approximate  solution  of  (2.13). 
By  assuming  a;^(y)  = for  constant  4>,  we  get 


a>  = 


1(0) 


» + 


•Okx, 


(3.1) 


2M0)|??,(k)|. 


where  ^k^  is  obtained  from  a formula  in 


“* " ' (3.2) 


and  - is  the  Strouhal  span  number. 


The  coefficient  of  the  lifting  force  of  the  wing  is  calculated 
by  formula 

^^  \ ya.  lOdld%  ^ (3.3) 
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Now  let  us  write  the  coefficient  for  the  lifting  force  of  a 
wing  in  the  form  of  two  components: 

— Cy,  4" 

where  C is  the  coefficient  of  the  lifting  force  determined  by 
ye 

vortex  movements:  C - the  coefficient  of  the  lifting  force 

yii 

determined  by  the  connected  mass  of  the  wing. 


Now  let  us  find  components  C,  and  C : 

^ ys  yw 


Here 


- ..  '.A  .-1 + 

+1 

C,.  = -%2ik  J VT^F,  (x)  dx. 


(3.^) 


(3.5) 


and  C(k)  is  the  Theodorsen  coefficient  [2]. 

In  calculations  the  theoretical  value  of  a = 2tT  (with  the 
effect  of  the  viscosity  of  the  fluid  considered)  can  be  assumed 
to  be  a =5.^5. 


For  a stationary  flow  = 0,  while 
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+ 


2»f-^xdxri  — j.,  1 


(3.6) 


The  first  term  in  (3.6)  gives  us  the  well  known  value  of  the 
wing  coefficient  in  a limited  flow;  the  second  term,  which  has  an 
order  of  1/X,  represents  the  refinement  which  the  present  method 
introduces  into  the  Prandtl  theory . For  a plane  lifting  surface 
the  second  term  in  (3-6)  is  equal  to  zero,  and  the  results  will 
be  in  complete  agreement  with  the  classical  theory  of  the  support- 
ing line  [3,  4]. 


Functions  (p^,  (p^,  p were  obtained  in  [7,  8,  11]  for  forward 
and  reciprocal  wing  oscillations. 
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